Note on log-periodic description of 2008 

financial crash 
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We analyze the financial crash in 2008 for different financial markets from the 
point of view of log-periodic function model. In particular, we consider Dow Jones 

index, DAX index and Hang Seng index. We shortly discuss the possible relation 

O 

CN ' of the theory of critical phenomena in physics to financial markets. 

a' 

o\ ; 1 Introduction 

^ . In recent years many papers have appeared which deal with applications of the 

S ' methods of statistical physics to the stock market description, in particular - to the 

qh' analysis and prediction of financial crashes. 

The interactions between physical particles are usually short-ranged. In spite 
CN ■ 
^ [ of that physical systems exhibit, under certain conditions, cooperative behaviour. 
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% 



The Onsager solution [T] of the twodimensional Ising model may be viewed as the 
first mathematically precise derivation of this fact from first principles. 



O I It is very tempting to pursue this idea when considering the stock markets 



dynamics; here also a cooperative phenomena appear, for example in the form of 
financial crashes which may arise if the players in bulk decide to sell their shares. 
Moreover, they are usually neither associated with specific new items nor caused 
by the fact that a large group of agents communicate with each other. On the 
contrary, the traders can imitate the opinions of rather restricted set of nearest 
neighbors ("short-range interactions"). The mechanism leading to the observed 
stock market behaviour has been proposed by Johansen et al [2]-^[S]- They assume 
that the traders influence each other only locally and either imitate the opinions 
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of their nearest neighbors or take decisions independently. In the former case, if 
the tendency for traders to imitate their neighbors increases up to a certain critical 
point, many traders may place the same order - sell at the same time thus causing 
a crash. In the latter the buyers and sellers disagree with each other and roughly 
balance out. These ideas can be formulated analytically providing quantitative rela- 
tionships between some relevant parameters like critical exponents and log-periodic 
frequencies (see Sec. 2). 

The resulting model has been applied to many financial crashes (including 
real-estate p], oil [7] bubbles as well as some Chinese stock market bubbles |!B])and 
still invokes controversies [H], [ID]- The main one concerns, of course, the possibility 
of prediction when the crash will actually occur. According to the model the crash 
time is the accumulation point of geometric series of local minima. However, the 
very recognition of proper minima to be taken into account is not always easy. Even 
if we find it possible the result may appear unsatisfactory. As an example consider 
the time series of Hang Seng index in period 2003-2008. It has four minima at the 
distances: t2 — ti = 512, t^ — t2 = 191, t^ — t-^ = 102. The value of scaling factor A, 
as defined by eq. (ITU]) below, reads A = 2.68 or A = 1.87, depending on the points 
selected. 

It is, therefore, an open question whether the analogy between the critical 
phenomena in physics and the behaviour of stock market is deep or rather superficial. 
We remain slightly sceptical. However, we believe that further analysis is necessary. 
In the present paper the analysis of financial crash in 2008 for different markets is 
presented. 

The paper is organized as follows. In Sec. 2 the theoretical model introduced by 
Johansen et al. is sketched. Sec. 3 contains empirical results. Last section is 
devoted to some conclusions. 
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Figure 1: Logarithm of Hang Seng index over the period 2003-2007. 



2 The JLS model 



The following theoretical model has been introduced by Johansen et al. [2] -^[5]. 
One considers the hazard rate which is the probability per unit time that the crash 
will happen in the next instant if it has not happened yet. Therefore, the hazard 
rate is defined as: 



i-QW 

where q (t) is the probability density function of the time of the crash, Q (t) - the 
cumulative distribution function. The dynamics of the asset price is given by: 

dp = fi (t) p (t) dt - up (t) dj. (2) 

The authors assume that, in the case of a crash, the price drops by a fixed percentage 
K. G (0, 1); j denote a jump process whose value is zero before the crash and one 
afterwards; /i (t) is chosen so that the price process satisfies the condition: 

W>t Et[p{t')]=pit). (3) 



This condition yields 



Inserting eq.(jl]) into eq. 



/i (t) = Kh it) 

we have 
pit) 



log 



Pito) 



K / h{t')dt' 
to 



(4) 



(5) 



before the crash. In order to be able to use eq. ([5]) the inforniation concerning the 
behaviour of hazard rate h (t) is needed. Johansen et al. [1] proposed the form of 
h (t) based on the above - mentioned analogy with statistical physics. Namely, they 
assumed that any agent i can be only in one of two possible states Si G { — 1,+1} 
(buy/sell). The state of i-th agent is determined by: 

Si = sign I K \^ Sj + aSi I (6) 

where the sign{-) function is equal to +1 (—1) for positive (negative) numbers. 
A^ (i) denotes the set of neighbours of i-th agent; only their decisions can influence 
him/her. K is positive constant, Ei are independently distributed according to the 
standard normal distribution. K is called coupling strength and control the ten- 
dency towards imitation; the tendency towards idiosyncratic behavior is governed 
by a. 

The model formulated in the above way resembles the Ising model solved (in twodi- 
mensional case) by Onsager. Thus one can expect that there exists a critical value 
of K, Kc, such that the system becomes very sensitive to small global perturba- 
tions (in the physical language, it gets magnetized after applying an arbitrarily 
small external magnetic field and exhibits large-scale fluctuations). It is assumed 
that the hazard rate which measures the reaction on small perturbations behaves 
as magnetic susceptibility i.e. 

h (t) ^ Bo {U -ty^ + Biitc-ty" COS [co log (t, -t)+ ij'] . (7) 

where the critical time tc is defined as the first time such that K (tc) = Kc 
Putting eq.([7]) into eq.([5]), one arrives at the following formula for the evolution of 
the price before the crash 

log [p (t)] ^A + B{tc-tf{l + Ccos [oo log (t, - t) + 0]} (8) 

where (p is another phase constant. 



3 Empirical Results 

3.1 Fitting the log-periodic function to empirical date 

The empirical search for log-periodic function describing the finacial market was 
carried out in many papers [lT]-[16]. In order to find seven parameters of log- 
periodic function we must minimize the following expression: 



tn t 



^ {yt - ytf = Y,{yt-A-B{t,- t)" (1 + Ccos (u;log (t, -t)+ ct>))f , (9) 
t=ti t=ti 

where yt is the logarithm of real index price, yt is the date point as predicted by 

the model, and t^ < t^. However, the minimalisation is not an easy task. The 

noisy data and a fitting function with large number of degrees of freedom cause 

many local minima of eq. ([9]) to exist where the minimalisation algorithm can get 

trapped. In some papers |1],[IS] the parameters A, B, C are reduced by requiring 

first that eq.([n]) has at the minimum the vanishing derivatives with respect to A, 

B, C. In this way one gets 3 linear equation for the latter. Solving them allows 

to express A, B and C as functions of remaining parameters. However, the form 

of these functions is very complicated and, moreover, one has still quite a lot of 

parameters controlling the fit. 

For this reason we look for the best fit in the following way. The log-periodic 

corrections to scaling imply the existence of hierarchy of characteristic scales in space 

or time. The linear scale is determined by three consecutive minima or maxima of 

log-periodic oscillations [T2].[T7] 

A = ^"+^~^" . (10) 

tn+2 ~ tn+1 

In eq. (|8]) the linear scale is contained in parameter u = 2t\ j In A. The authors of 
papers [II]-[13] postulated that A is the universal parameter approximately equal 
2. However, eq. ( ITOl) shows that the A parameter may in principle depend on the 
specific geometry and structure of the system. 

Assuming that the spacing beetwen successive values of t„ approaches zero as 
n becomes large and t^ converges to t^ (more precisely, assuming eq. (ITU]) to hold 



with A > 1), we can show that: 

2tn+l — tfi — tn+2 

Thus we determine A and tc from three successive observed values of t„. There 
remain five parameters which are then determined with the help of nonlinear re- 
gression (Nonlinear Regression packet of Mathematica). This preliminary fit allows 
us to put some constrains on parameters we are looking for. Using these constraints 
we minimize the right hand side of eq. iQ with respect to all seven parameters; due 
to the constraints obtained in the first step it is now unlikely to get false minima. 

3.2 Dow Jones Industrial Average Index 

We analyze the bubble data for the period from 02-Feb-2003 to 31-Oct-2007. One 
must keep in mind that the system is scale invariant only nearby the critical point. 
For this reason we carry out the procedure of fitting the log-periodic function to 
real data only for data between ten and twenty months before the crash. From the 
period 18.01.2006 - 29.08.2007 three successive minima have been selected in order to 
calulate A and tc- From equations ( ITOj) and (TTT]) we find A = 1.4, and tc = 18.09.2008 
(which corresponds to the point 672). We looked for the parameters of the best fit 
of the log-periodic function: 

log [p (t)] ^A + B (672 - t)" {1 + C cos [18.66 * log (672 - t) + 0]} (12) 

to real data. The values of these parameters are presented in Table [H The value of 
estimated variance, defined as MSE = {J2t=t (Vt ~ Vt) ) Z*^/' ^^ equal to 3.68 x 10~^; 
here df = 402 is the number of degrees of freedom. 

The parameters estimated from nonlinear regression do not allow to determine 
very restrictive constraints for all parametrs of eq. (IHl)- We assume that 15 < 
w < 20, 5 < 0, < < 27r, 0.1 < a < 1, 06.06.2008 < tc < 28.10.2008. The 
ranges of acceptable values for parameters are greater than the confidence intervals 
(Table 1) because we haven't estimated the values of A and tc- With the above 
constraints the function i^ was minimalized. We obtained the following values of 
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Table 1: The values of estimated parameters of log-periodic function for Dow Jones 
index from 18.01.2006 - 29.08.2007. 





Estimate 


Asymptotic SE 


Confidence Interval 


A 


8.68961 


0.0659957 


{8.55987,8.81935} 


B 


-0.00405699 


0.00410804 


{-0.0121329,0.00401894} 


C 


0.0565023 


0.0114346 


{0.0340232,0.0789814} 


a 


0.737919 


0.13599 


{0.470578,1.00526} 





2.45213 


0.0632643 


{2.32776,2.5765} 



the parameters: A -^ 8.8106, B -^ -0.0165957, C -^ -0.0444881, a -^ 0.554188, 
tc —7- 672.319, (/?—)■ 0, a; — 7- 19.5637, and the estimated variance is 3.51 x 10""^. 
The determination coefficient equals 93% and is a slightly greater than for the case 
when A and t^ were calculated from eq. fITII]) and flTTl) . It is very interesting that 
the estimated date of crash is the same as the one which has been calulated from 
eq. (TTTI) . Figure |3] shows that the estimated log-periodic function describes very 
well the price dynamics in the whole period under consideration, despite of the fact 
that it was found by considering a short period (Figure [2]). 

Log(index) 
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Figure 2: Logarithm of Dow Jones index over the period 18.01.2006 - 29.08.2007 
versus its corresponding log-periodic representation. 
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Figure 3: Logarithm of Dow Jones index over the period 12.11.2002 - 26.10.2009 
versus its corresponding log-periodic representation. 

Table 2: The values of estimated parameters of log-periodic function for DAX index 
from 25.05.2006 - 27.08.2007. 





Estimate 


Asymptotic SE 


Confidence Interval 


A 


9.4309 


0.0988434 


{9.23643,9.62538} 


B 


-0.0158149 


0.0118548 


{-0.0391392,0.00750939} 


C 


0.0440315 


0.00673482 


{0.0307807,0.0572823} 


a 


0.632698 


0.101428 


{0.433139,0.832257} 


^ 


0.245255 


0.074815 


{0.0980567,0.392454} 



3.3 DAX Index 

DAX index measures the development of the 30 largest and best-performing compa- 
nies on the German equities market and represents around 80% of the market capital 
authorized in Germany. From the data concerning the period 25.05.2006-27.08.2007 
one finds using eqs. ( ITOj) . (TTT]) A = 1.55 and t^ = 26.05.2008. Other estimated pa- 
rameters are shown in Table 121 The estimated variance is equal 6.44 x 10~^. On 
the basis of these values, we assume the following constraints: 03.03.2008 < tc < 
30.12.2008, 10 < w < 18, 5 < 0, 0.1 < a < 1. The minimalization of eq. (Q 
gave the following results: A -^ 9.52346, B -^ -0.0130455, C -^ -0.0390518, 
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a -^ 0.669289, r -^ 570.63 ^ 21.08.2008 (actually, the crash occured in October; 
however, the DAX index value fell by above 17% already in the period 08.-10.2008), 
Lf — > —2.33041, u — )■ 18. The estimated variance is equal to 6.026 x 10"^ which 
means that the regression curve approximates the real data in 96%. In contrast 
to Dow Jones the estimated log-periodic function of DAX does not describe very 
well the price dynamics in the whole period (Figure S]). It reflects, however, the 

Log(index) 
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Figure 4: Logarithm of DAX index over the period 25.05.2006 - 27.08.2007 versus 
its corresponding log-periodic representation. 
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Figure 5: Logarithm of DAX index over the period 29.09.2003 - 24.09.2009 versus 
its corresponding log-periodic representation. 



decrease of price in August 2004 and October 2006. 



3.4 Hang Seng Index 

The Hang Seng Index is used to record and monitor daily changes of the largest 
companies of the Hong Kong stock market and is the main indicator of the overall 
market performance in Hong Kong. 

In this case we decided to check whether it is possible to get a reasonable description 
by taking into account the data from a longer period, not from the one directly 
preceding the crash. In statistical physics the universal behaviour is expected in 
the neighbourhood of critical point while far from it the properties of the system 
may be very specific and depend on particular dynamics. However, it appears here 
(see below) that the fit with the help of log-periodic function is still satisfactory. 

We take four successive minima from which, from eq. flTOj) and f fTTj) . we get the 
average value of parameter A ~ 2.3 and t^ = 1087 — ?■ 18.02.2008. Using them, we 
put the constraints: 990 < tc < 1200, 7 < a; < 11, < < 27r, fi < 0, 0.1 < a < 1. 
The parameters minimalizing eq. (|9]) read: A — )■ 11.3902, B — )■ —0.550977, C -^ 
-0.00857998, a -^ 0.184217, t^ -^ 1085.13, ^ 0, w ^ 8.90061. The regression 
function approximates the real data points as much as in 96% (Figure E]); the system 
had 972 degrees of freedom. 




date 
07.2004 05.2005 02.2006 12.2006 10.2007 08.2008 

Figure 6: Logarithm of Hang Seng index over the period 29.09.2003 - 09.09.2009 
versus its corresponding log-periodic representation. 
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3.5 The next terms of log-periodic function 

So far we were considering the log-periodic function only keeping the first term of 
Fourier expansion in eq. ([8]). Now, we analyze the price dynamics of Dow Jones 
index taking into account both the first and second term of Fourier expansion. 
We choose the same time period as subsection 3.2. The log-periodic function is 
described by: 

log [p it)] ^ A+B {tc -t)"{l + C cos [u log (tc -t) + (j)]+D cos [2u log (t^ - t) + ^]} 

(13) 
With the use of the same constraints on parameters as in subsection 3.2, we have: 
A -^ 8.93833, B -^ -0.0423764, C -^ -0.0367758, D -^ -0.0173993, a -^ 
0.439734,te -^ 667.641 ^ 12.09.2008, i) -^ -8.64053, ^ 0, w ^ 19.6062. The 
regression curve approximates the real data in 94,2% which, as expected, is more 
than in the case when only the first term is taken into account. The approximated 
function describes the price dynamics of bubbles better but the improvement is not 
dramatic (Figure [7]). We could expect that the next terms of Fourier expansion of 
log-periodic function will increase the determination coefficient but not significantly. 

Log(index) 




date 



Figure 7: Logarithm of Dow Jones index over the period 12.11.2002 - 26.10.2009 
versus its corresponding log-periodic (first and second terms of Fourier expansion) 
representation. 
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4 Conclusions 

We analyzed the financial crash in October 2008 on the basis of indices for different 
markets. Only in the case of Hang Seng index we received the value of critical 
time compatible with the definition of critical point as the point where the central 
trend of price dynamics changes its direction [T3]. In other cases it was rather the 
moment when the index price drops sharply in short period. This finding supports 
the scepticism concerning the possibility of exact prediction of crash date on the 
basis of log-periodic function approach expressed already in Refs. 0, [lU], [IB]- 

The same financial crashes were discussed in the papers [TU] and [13]. Their 
authors obtained A = 2 and A = 3, respectively. The differences seem to result 
from the way the minima of time series were selected. Therefore, we think that 
the predictive power of the approach based on log-periodic functions is doubtful. 
However, the aposteriori description of the price history is quite satisfactory as it is 
seen from the large value of determination coefficient. 

It is very tempting to compare the djTiamics of market indices with the crit- 
ical phenomena in condensed matter physics. The latter offer the well understood 
example of the situation when local interactions lead to large-scale cooperative 
behaviour; the main idea is the emergence of scaling invariance, the main tool - 
the renormalization group transformations which provide a natural framework for 
studying the universality phenomena. However, it is rather unlikely that the expla- 
nation offered by the theory of critical phenomena is sufficiently general to cover all 
cases of emergence of cooperative behaviour in complex dynamical systems. 

The analogy with statistical mechanics can be also questioned due to the fact 
that, contrary to the case of condensed matter systems, the description based of 
scaling seems to work reasonably also quite far from critical points. As it had been 
shown above, in the case of Hang Seng index the determination coefficient was equal 
96%. Guided by this value we have verified also the quality of approximation for the 
data from longer period 2003-2008 for remaining indices. The relevant coefficients 
for Dow Jones and DAX were equal 97.47% and 97%, respectively. Therefore, the 
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distance from the criticality region does not seem to be a really relevant parameter, 
contrary to what is expected in physics. 
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